Quantum Hall effect beyond the linear response approximation 
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The problem of Bloch electrons in two dimensions subject to magnetic and intense electric fields is 
investigated, the quantum Hall conductance is calculated beyond the linear response approximation. 
Magnetic translations, electric evolution and energy translation operators are used to specify the 
solutions of the Schrodinger equation. For rational values of the magnetic flux quanta per unit 
cell and commensurate orientations of the electric field relative to the original lattice, an extended 
superlattice is defined and a complete set of mutually commuting space-time symmetry operators are 
obtained. Dynamics of the system is governed by a finite difference equation that exactly includes the 
effects of: an arbitrary periodic potential, an electric field orientated in a commensurable direction 
of the lattice, and coupling between Landau levels. A weak periodic potential broadens each Landau 
level in a series of minibands, separated by the corresponding minigaps; additionally the effect of 
the electric field in the energy spectrum is to superimpose equally spaced discrete levels, in this 
"magnetic Stark ladder" the energy separation is an integer multiple of hE/aB, with a the lattice 
parameter. A closed expression for the Hall conductance, valid to all orders in E is obtained, the 
leading order term reduces to the result of Thouless et al, in which tj^' is quantized in units of 
e'^/h. The first order corrections exactly cancels for any miniband that is completely filled. Second 
order corrections for the miniband conductance are explicitly calculated as cr^' oc e^/USB, with Uo 

(2) 

the strength of the periodic potential. However the use of a sum rule shows that cancels when 
a Landau band is fully occupied. 



I. INTRODUCTION 

The quantization of the Hall conductance (cfh) in 
integer-i and fractional values'^ is an astonishing and un- 
expected result, because in a real sample there are many 
effects that would be expected to flaw this possibility. 
In the integer quantum Hall effect (IQHE), the conduc- 
tance is an integer multiple of /h; simultaneously the 
longitudinal magnetoresistance vanishes; both results are 
observed within an experimental uncertainty of less than 
3 X 10~^. For the IQHE the exact quantization of the 
Hall conductivity was explained based on a topological 
reasoning by Thouless et almi^; an in units of /h is a 
topological invariant, the so-called first Chern numbesSii, 
which can only take integer values, ah is an integer mul- 
tiple of /h if the Fermi level lies in a gap between Lan- 
dau levels. The addition of a periodic potential broadens 
each Landau level in a series of minibands, separated 
by the corresponding minigaps. Contrary to the obvi- 
ous suggestion that each subband carries a fraction of 
e^//i, Thouless et a/.demonstrated that each filled mini- 
band also contributes with an integer multiple of e^ /h to 
gh- The Hall conductance varies in a non-monotonous 
sequence as the Fermi level sweeps through contiguous 
minigaps. Clear experimental evidence for the internal 
structure of the Hofstadter butterfly spectrum has only 
been found very recently in the measurement of uh for 
lateral superlattices?. The work of Thouless et aLmakes 
use of the Kubo linear response theory, it is of our inter- 
est to consider the structure of the Hall current when a 
finite electric field [E) is applied, in particular we shall 



be able to compute higher order corrections on E to the 
Hall conductivity. 

The problem of electrons moving under the simultane- 
ous influence of a periodic potential and a magnetic field 
has been discussed by many author3Sii2iiiii2iiii^iiSii&; 
the spectrum displays an amazing complexity including 
various kinds of scaling and a Cantor set structureii. 
Some of these results had been used in order to disclose 
the topological structure of the Hall conductance within 
the linear response theory. The addition of the electric 
field leads to the possibility of testing the exactness of 
the quantization beyond the linear approximation. Fur- 
thermore the problem by itself possesses a rich physical 
structure that makes its analysis worthwhile. 

We consider the problem of an electron moving in a 
two-dimensional lattice in the presence of applied mag- 
netic and electric fields. We refer to this as the electric- 
magnetic Bloch problem (EMB). The corresponding 
magnetic Bloch system (MB) has a long and rich history. 
An early important contribution was made by Peierla^ 
who suggested the substitution of the Bloch index k by 
the operator (p — eA) in the B = Q dispersion relation 
e(fc), which is then treated as an effective single-band 
Hamiltonean. The symmetries of the MB problem were 
analyzed by Zakiiii^ti^, who worked out the representa- 
tion theory of the of magnetic-translations group. The 
renowned Harper equation was derived assuming a tight- 
binding approximation in which the magnetic field acts 
as a perturbation that splits the Bloch bands^". Rauh 
derived a dual Harper equation in the opposite limit of 
intense magnetic fieldii, here the periodic potential acts 
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as a perturbation that broadens the Landau levels, in this 
case the Harper equation takes the form 



II. THE ELECTRIC-MAGNETIC BLOCK 
PROBLEM 



+ 2 cos (27rf7m + k) c„i + c„i+i 



(1) 



The model 



where a = \/<j) is the inverse of magnetic flux </> in a 
cell in units of h/e. The studies of the butterfly spec- 
trum by Hofstadter and others have since created an un- 
ceasing interest in the problem because of the beautiful 
self-similar structure of the butterfly spectrumii. Re- 
markably, an experimental realization of the Hofstadter 
butterfly was achieved considering the transmission of 
microwaveSii^ and acoustic waves^S through an array of 
scatterers. 

The symmetries of the EMB problem were analyzed 
by Ashby and Miller, who constructed the group of the 
electric-magnetic translation operators, and worked out 
their irreducible representationE^*^. The properties of the 
electric-magnetic operators were utilized in order to de- 
rive a finite difference equation that governs the dynamics 
of the EMB problem when the coupling between Lan- 
dau levels can be neglectedli. In this paper we shall 
derive the equation that applies under most general con- 
ditions. Magnetic translations, electric evolution and en- 
ergy translation operators are used to specify the solu- 
tions of the Schrodinger equation, commensurability con- 
ditions must be implemented in order to obtain a set of 
mutually commuting space-time symmetry operators. In 
addition to the broadening of the Landau levels induced 
by the periodic potential, the effect of the electric field 
in the energy spectrum is to superimpose equally spaced 
discrete levels; in this "magnetic Stark ladder" the en- 
ergy separation is an integer multiple of hE/aB, with a 
the lattice parameter. A closed expression for the Hall 
conductance, valid to all order in E is obtained. The 
leading order contribution is quantized in units of e? /h. 
First and second order correction are explicitly computed 
and shown to take the form: cr^-' — and cr^-* cx /U^B 
respectively, for any miniband that is completely filled. 

The paper is organized as follows. In the next section 
we present the model that describes the EMB problem 
and construct its symmetry operators. In Section IIIII 
we describe the commensurability conditions required to 
have simultaneously commuting operators, they are ex- 
ploited in order to construct an appropriated wave func- 
tion basis. We derive an effective equation (|30|) . in which 
the "evolution" is determined by a differential equation 
with respect to the longitudinal pseudomomentum, this 
equation becomes essential in order to obtain a closed ex- 
pression for the Hall conductance (|74|) . The derivation of 
the finite difference equation that governs the dynamics 
of the system is presented in Section Hvl Results for the 
energy spectrum are presented, and are also discussed 
from the perspective of the adiabatic approximation. A 
closed expression for the Hall conductance, valid to all 
order in E is obtained in^ the leading order term cr^' 

and perturbative corrections cr^'' and cr^^ are calculated. 
The last section contains a summary of our main results. 



Let us consider the motion of an electron in a two- 
dimensional periodic potential, subject to a uniform mag- 
netic field B perpendicular to the plane and to a con- 
stant electric field E^ lying on the plane according to 
E = E {cos 9 , sin 9) with 6 the angle between E and the 
lattice a:— axis. The dynamics of the electron is governed 
by a time-dependent Schrodinger equation that for con- 
venience is written as 



1 



2m* 



Hi 



n^) + y - Ho 



1^) = 0, (2) 



here m* is the effective electron mass, H^ = + e^^, 
with — [ihd/dt^—ihV). Except for the final re- 
sults for the conductivity, throughout the paper energy 
and lengths are measured in units of hujc = and 

If) — \J respectively, where u>c is the cyclotron fre- 
quency and i{) the magnetic length. So, unless specified, 
we set h = e = m* = 1; although B can also be omitted 
from the expressions, we find convenient to explicitly dis- 
play it. Covariant notation will be used to simplify the 
expressions, e.g. = {t,x) — {t,x,y). 

Equation (|2Jl can be considered as an eigenvalue equa- 
tion for the operator S with eigenvalue 0. The gauge 
potential is written in an arbitrary gauge, it includes 
two gauge parameters a, and (3, the final physical results 
should be of course, independent of the gauge. Hence the 
components of the gauge potentials are written as 



/3+-]E.,t 



f3- 



E„t - 



a - ^] By, 



a + - \ Bx . 



(3) 



The symmetrical gauge is recovered for a = 0, whereas 
the Landau gauge corresponds to the selection a = 1/2. 
If /3 = 1/2 all the electric field contribution appears in the 
vector potential, instead if /3 ~ —1/2 it lies in the scalar 
potential and the Schrodinger equation becomes time in- 
dependent, however even in this case a time dependence 
will slip into the problem through the symmetry opera- 
tors. A general two-dimensional periodic potential can 
be represented in terms of its Fourier decomposition 

/ 27rrx 2Tr sy\ 
V [Xi y) ^ '"rs exp I h I . (4) 

r,s ^ ' 

For specific numerical results we shall use the potential 

'^Uacos^^^ 



V{x,y)^Uf) 



cos 



\ a J 



(5) 
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A is a parameter that can be varied in order to have an 
anisotropic lattice; A = 1 corresponds to the isotropic 
limit. 



III. ELECTRIC-MAGNETIC BLOCK 
FUNCTIONS 

A. Commensurability conditions 



B. Electric evolution and magnetic translations 

Let {t, x) {t + T,x + R) be a miiform translation 
in space and time, where r is an arbitrary time and R 
is a lattice vector. The classical equations of motion 
remain invariant under these transformations; whereas 
the Schrodinger equation does not, the reason being the 
space and time dependence of the gauge potentials. Nev- 
ertheless, quantum dynamics of the system remain invari- 
ant under the combined action of space-time translations 
and gauge transformations. The electric and magnetic 
translation operators are defined as 

To(t) = exp {-itOo) , Tj{a) = exp {laOj) , (6) 

with j = x,y and the electric-magnetic symmetry gener- 
ators are written as covariant derivatives = p^^ -\- , 
with the components of the dual gauge potentials 
given by 



Ao = An+x- E , 

K = +By + E^t , 



Ky = Ay — BX 



Eyt. 



(7) 



It is straightforward to prove that the operators in © 
are indeed symmetries of the Schrodinger equation; they 
commute with the operator S in Eq. Q. Similar ex- 
pressions for the electric-magnetic operators were given 
by Ashby and Miller^°, however their definition included 
simultaneous space and time translations; we deemed it 
more convenient to separate the effect of the time evo- 
lution generated by the Tq to that of the space transla- 
tions generated by Tj. The following commutators can 
be worked out 



[no,n,] 
[n^,a] = 0. 



lE, 



[ni,n2] 



-iB, 
iB, 



(8) 



Notice that the electric-magnetic generators have 
been defined in such a form that they commute with all 
the velocity operators 11^. The Schrodinger equation and 
the symmetry operators are expressed in terms of covari- 
ant derivatives and , respectively. A dual situation 
in which the roles of 11^ and are interchanged can be 
considered. The dual problem corresponds to a simulta- 
neous reverse in the directions of B and E. We notice 
that the commutators in the second line of Eq. ((HJ are 
part of the Lie algebra of the EM-Galilean two dimen- 
sional group22i2^. This group is obtained when the usual 
rotation and boost operators of the planar-Galilean group 
are replaced by their electric-magnetic generalization in 
which the operators are enlarged by the effect of a gauge 
transformation. 



In order to construct a complete base that expands 
the wave function we require the symmetry operators to 
commute with each other. However we have 



J- J J J- II. 



(9) 



A set of simultaneously commuting symmetry operator 
can be found if appropriated commensurate conditions 
are imposed, we follow a three step method to find them: 



1, 



First we consider a frame rotated at angle 9, with 
axis along the longitudinal and transverse direc- 
tion relative to the electric field. An orthonormal 
basis for this frame is given by = {cos6,sin9), 
bt = (—sinO, cos9) and 63 = x e^. The electric 
field is parallel to and the magnetic field points 
along 63. We assume a particular orientation of 
the electric field, for which the following condition 
holds 



tan 9 



Ey 

Err 



TO2 



(10) 



where mi and m2 are relatively prime integers. 
This condition insures that spatial periodicity is 
also found both along the transverse and the lon- 
gitudinal directions. Hence, we define a rotated 
lattice spanned by the longitudinal 6^ = be^ 
and transverse bx — hex vectors, where b = 
a-y/ mf + m\. The spatial components of the sym- 
metry generator C are projected along the longi- 
tudinal and transverse directions: = e^, • C 
and Ot = ex ■ O. It is readily verified that 
[Oo,Ot] =0. 

2. For the rotated lattice, we regard the number of 
flux quanta per unit cell to be a rational number 
p/q, that is 



1 
a 



B9 
27r 



P 

q 



(11) 



We can then define a extended superlattice. A rect- 
angle made of q adjacent lattice cells of side b con- 
tains an integer number of flux quanta. The basis 
vectors of the superlattice are chosen as qb^ and 6t- 
Under these conditions the longitudinal and trans- 
verse magnetic translations {qb) = exp{iqbOL) 
and Tt (6) = exp [ibOr) define commuting symme- 
tries under displacements qb^ and by- Henceforth 
we shall either use the subindex {L,T) or (1,2) to 
label the longitudinal and transverse directions. 

3. We observe that Tq and TL{qb) commute with Tt- 
Yet they fail to commute with each other: 



To (r) Tl {qb) = e-'^^-^Ti {qb) Tq (r) 



(12) 
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However the operators Tq and (qb) will commute 
with one another by restricting time, in the evolu- 
tion operator, to discrete values with period 

r = nro, n G Z, tq = ^ - [ — ] , (13) 

qbE p \VD J 

where the drift velocity is vu = E/B and we uti- 
lized Eq. (fTT|l to write the second equality. b/v£, 
is the period of time it takes an electron with drift 
velocity vd to travel between lattice points. Eq. 
H13|l can be interpreted as the condition that the 
ratio of two energy scales is an integer: as we 
shall see (abE) is the magnetic-Stark ladder spac- 
ing, whereas {2'iTVD/h) is the quasienergy Brillouin 
width; hence Eq. (|13|) represents the ratio of these 
two quantities. 

We henceforth consider that the three conditions pO|l . 
and (|13|l hold simultaneously. In this case the three 
EM operators: the electric evolution Tq = To(to), and 
the magnetic translations 7^ = Ti(g6), and 7^ = Txip) 
form a set of mutually commuting symmetry operators. 
In addition to the symmetry operators it is convenient to 
define the energy translation operatOB24 



Tg — exp —t — t 

To 



(14) 



that produces a finite translation in energy by 2tt/tq = 
qbE. Te commutes with the three symmetry operators 
but not with S. Its eigcnfunctions 



define a quasitime d modulo tq. 



(15) 



B. Wave function and generalized Bloch conditions 

Having defined Tq , Tj, and that commute with each 
other and also with S*, it is possible to seek for solutions 
of the Schrodinger equation labeled by the quasienergy 
(£) and the longitudinal (ki) and transverse (^2) quasi- 
momentum according to 

ro|£,fci,fc2> =e-'^"^|£,fci,fc2), 

rL|£,fci,fc2> =e''^'?^|£,fci,fe), (16) 

Tt If, fci, fc2) = 6^'=^'' |f,A:i,fc2). 

The magnetic Brillouin zone (MBZ) is defined by fci G 
[0, 27r/q6] and fc2 € [0, 27r/6]. Similarly the quasienergy 
£ is defined modulo 27r/ro = qbE. If a restricted energy 
scheme is selected for the energy, the first energy Bril- 
louin region is defined by the condition £ £ [0,2tt/to]. 
We shall find convenient to perform a canonical transfor- 
mation to new variables according to 



Qo=t, 

Ql =H2 



Pn = a - Ie 



Et, 



Pi =111, 

P2=02 



(17) 



E, 



that satisfy the commutation rules [Qfj.,Pu] = 
—ig^v, g^iv = -D*aff(— 1, 1, !)• Apphed to Eq. Q the 
transformation yields for the Schrodinger equation 



"1 



H 



Pt + Qi)+V{x,y)-EP2 



(18) 



In this equation (x, y) must be expressed in terms of the 
new variables: 

^=^(Qi-p.)-^(g.-po, 

a b 

^ (Qi-P2) + -^ (Q2-P1) . (19) 

On the other hand, the EM-symmetry operators and the 
energy translation operator take the form: 



■J- = ^-'iroPo 



io = e 

n- ibP2 

It = e ^ 



7—1 — 
E = e ^0 



Notice that these operators do not depend on the vari- 
ables Pi,Qi. Thus, it is natural to split the phase 
space {Qf,,Pf_c) in the ((5i,Pi) and the {Qo, Po;Q2, P2) 
variables. For the first set a harmonic oscillator base 
is used. Whereas for the subspace generated by the 
variables (Qo, Po! Q2, P2), we observe that the opera- 
tors {Tq, Tl, Trj^, T^) with all possible integer values of 
{i,i,k,l) form a complete set of operators. The demon- 
stration follows similar steps as those presented by Zak 
in references^. Hence a complete set of functions, for 
the subspace (Q07 Po; Q2, P2), is provided by the eigcn- 
functions of the operators {T^, Tl, Tr^, T^). As starting 
point, we select a base of eigenvalues of the following 
operators 



A^A £, k2) = n IfJ., £, k2) , 

Po \fJ;£,k2) = £ \fJ.,£,k2) , 

P2 \n,£,k2) = k2 \n,£,k2) , 



(21) 



where A and A'^ are the lowering and raising operators 
of the /i-Landau level: 
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At = — (Pi+iQi). (22) 



It is straightforward to verify that these states fulfill the 
required eigenfunction condition (|16|l for the Tq and Tx 
operators. On the other hand Tl induces a shift in 
the £ and fc2 labels, while Te induces a shift qbE in 
the £ label. The previous considerations suggest that 
a state \'d,£,ki,k2) that is a simultaneous eigenfunc- 
tion of the four operators Te, Tq, Tl and Tt with the 
required eigenvalues (Eqs. El £^nd ITCI) . can be con- 
structed as a linear superposition of states of the form 
[Tl]'' \TlTe'^Y^ \ii,£,k2). We write down the state, and 
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verify their correctness: 

\d,E,ki,k2) = 



(23) 



It is easy to check that this function satisfies the three 
eigenvalue equations for the symmetry operators in (|16|l . 
Additionally it can be verified that the eigenvalue condi- 
tion H15|l for the energy translation operator is also ful- 
filled by imposing the periodicity condition c^^p = c^. 
Every state in H23(l yields a set of different eigenvalues, a 
condition that follows from the fact that the selected set 
of operators is complete, hence the base of eigenvalues in 
(|23|l is complete and orthonormal. 

Up to this point we have constructed a base for the 
set of four EM translation operators, however the oper- 
ator Te is not a symmetry of the problem, consequently 
the solution of the Schrodinger is constructed as a su- 
perposition of all the states characterized by 1!); the state 
becomes 



I 



where 



(25) 



Equation (|24|l represents the correct state that satisfies 
the eigenvalue equations It is convenient to recast 

it in a compact form as 



\£,k)^W{ki)\£,k2), 



(26) 



where |f , k) — |£, /ci, ^2), and the operator W (fci) is de- 
fined as 



tqbl(EQo+Q2-ki) 



(27) 



whereas the ket \£,k2) is given by 

|£, k2) = e^^^™('3^-'=^)6^^ Im, £. k2) . (28) 

The previous expression stresses the fact that the quan- 
tity e^^'^bmki f^/j^ does not depends on ki, this will be 
demonstrated below Eq. H45I) . These results will prove 
to be very useful to analyze the EM problem, in par- 
ticular they allow us to obtain an effective Schrodinger 
equation where the "dynamics" is governed by the deriva- 
tive with respect to the longitudinal pseudomomentum. 
Utilizing Eqs. (|26|) y (|27l) it is straightforward to prove 
the following relation 

Pa\£,ki,k2) ^[Po,W]\£,k2)+WPo\£,k2) = 

{£ -iE-Vk)\£,ki,k2). (29) 



Thus the Schrodinger (|18|l can be recast as 

{£-tE- Vfc) \£, ki,k2) = H \£, fci, fe) . (30) 

As it will be lately discussed, this form of the Schrodinger 
equation becomes essential in the discussion of the quan- 
tized Hall current. 

The wave function takes a simply form if we adopt 
the (Po,Pi,P2) representation. In this case (p) = 
{Po,Pi,P2\£,ki,k2) yields 



J2 &^:,0^(Pi)e'^2''/''^('"+f')'=^ 



2tt 



5{Po-£- IqbE) S \^P2 -k2 + —im+pl)j , 

(31) 

where (Pi) is the harmonic oscillator function in the 
Pi representation 



(Pi) = (Pi|^) 



1 



(Pi) , (32) 



and Hfj^ (Pi) is the Hermite polynomial. 

On the other hand, it is sometimes useful to restore 
the space-time representation of the wave function. This 
problem provides a good example of the use of canonical 
transformations in quantum mechanics. Identifying the 
matrices that relate the (x^,p^) variables to (Q^, P^), it 
is possible to apply the method of reference^^ in order 
to obtain the desired transformation. The final result for 
the wave function in the space-time representation yields 



it,x) 



-iSt 



(33) 



where the modulation function u{t, x) is given by 



Uk (t, x) 



^-txi[{a~l/2}x2+ki]^t{P+l/2}Etxi ^ 
-iqbElt iab{ki —X2){Tn-\-pl) 



711+ pi 
b 

(34) 



here (f)^ is the same harmonic oscillator function Eq. H32|l . 
but now evaluated in the space representation. 

Solution (|33|) includes a superposition of Landau-type 
solutions originated beneath the spatial and time peri- 
odicity. The spatial periodicity is simply related to the 
external potential V . Whereas time periodicity arises 
from the conditions imposed to the symmetry operators 
in order to produce commuting symmetries; as discussed 
in Eq. l(T^ . the period is given by the time that takes 
an electron to drift between contiguous lattice points. 
Notice that Eq. H33|) follows from the Bloch and Floquet 
theorem. However in the electric-magnetic case the mod- 
ulation functions u (t, x) are not strictly periodic, instead 
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they satisfy the generahzed Bloch conditions 

U{t + To, Xl,X2) = e^'^"^" u{t, Xl,X2), 

u{t, xi + qb, X2) — e"*'^"^^ u{t, xi, X2), (35) 
u{t, xi,X2 +6) = e^^''^^ u{t, xi,X2), 

the phases are determined by the dual-gauge potential 
that appears in the symmetry operators Eq. 0. It is 
straightforward to verify that in Eq. H34|l satisfies 
these conditions. The correct normalization conditions 
for the modulation function are obtained as follows 



(2^)^ 
g&2 



cPx u* (t, x) u {t, x) = 1 , 



(36) 



MUC 



where MUC represents the magnetic cell defined by: xi S 
[0,qb] and X2 € [0,6]. 

The function u satisfies conditions similar to those in 
Eq. (I23J| but in the MBZ. The MBZ is actually a torus 
T^, so the edges (fci = 0, fc2) and (fci = 27r/g6, fc2) must 
be identified as the same set of points, the wave function 
can differ at most by a total phase factor (similarly for 
the edges (fci,fc2 — 0) and (fci,A:2 = 27r/6)): 



u{ki + 27r/(j6, k2) = e^^^uiki, k2), 
u{ki, k2 + 27r/6) = e'f^u{ki, fc2). 



(37) 



the functions /i (k) and /2 (fe) are related with the Hall 
conductance, however instead of Eq. (|35|l . these func- 
tions are not analytically known and must be computed 
numerically. 



IV. HARPER GENERALIZED EQUATION 

A. Finite difference equation 

The coefficient in Eq. (|24|) satisfies a recurrence re- 
lation that is obtained when this base is used to calculate 
the matrix elements of the Schrodinger equation (|18|l . 
First, let us consider the contribution arising from the 
periodic potential in Eq. Q). The x and y coordinates 
are written in terms of the new variables (Qi, Pi,Q2, P2) 
by means of Eq. Ijiyi) . producing a term of the form 
eyij>[(i2'Krmi{Qi — P2) /b], and another contribution in 
which Qi and P2 are replaced by Q2 and Pi. Once that 
Qi and Pi are replaced by the raising and lowering op- 
erators in Eq. (|22|l . one is lead to evaluate the matrix 
elements of the operator D = exp (zA^ — z*A) that gen- 
erates coherent Landau states. A calculation yields 

D^^" [z) = {v I exp {zA^ - z*A) \ fi) = 

yi^r" (1^1') . M > 
"S^r^ (1^0 ' < 



(38) 



where are the generalized Laguerre polynomial. The 
evaluation of the terms that include the P2 operator is di- 
rect, because the base is a eigenvalue of this operator, on 



the other hand the operator exp [{i2TTrmiQ2)/b] acts as a 
translation operator that produces a shift on the index m 
of the blfli coefficient. Taking into account these results, 
it is possible to demonstrate after a lengthly calculation 
that the Schrodinger equation (|I8I) becomes 



t=N 

E 

t=-N 



(fci , fca) kn+t - (f + Ek2 + abEm) L . (39) 



Besides its dependence on the index m, 6^ is a vector 
with L components and (fci, fc2) is a L x L matrix in 
the Landau space, according to 



r.s 



(40) 



In the previous expressions L is the highest Landau level 
included in the calculations, N — max{{r, s) (mi -I- 7712)} 
is related to the largest harmonic in the Fourier decom- 
position of the periodic potential |0J, and B'^{r,s) is 
given by 



(foo + A* + i) 



r, s = 0, 



(41) 

the following definitions where introduced in the previous 
equation 



(m2 — imi) [ir + s) , 
(rm2 — STOi) (rmi -\- sm2) , (42) 
{rmi + STO2) . 





V2tt 

b 




27r2 


Krs = 


^' 






Mrs = 









We recall that the indexes (r, I) refer to the Fourier 
expansion of the periodic potential Eq. whereas 
(mi, m2) correspond to the integers that determined the 
electric field orientation in Eq. H10|l. 

Equation H39|) describes the system in terms of a re- 
currence relation on a two dimensional (/i, m) basis. The 
coupling in m ranges from m — N to m + N. We fol- 
low the method of referenceSl to recast the recurrence 
relation on m with N nearest-neighbor coupling into a 
tridiagonal vector recurrence relation. This is enforced 
by defining a N component vector 



JNm, OjVm-l-1, 



-'Ar(m+1)-1 



, (43) 



7 



and matrices and Qm with elements 







-2 



V 



{ki,k2) = 



A 



-1 

Nm+l 



A 



\0 



y2~N 



Nm 



.N-2 



Using the fact that Q^^ = {Q^Y , Eq. can be reor- 
ganized as a tridiagonal recurrence relation 



Ek2 



f crbEm) In 



abEDN] 5,„ . (45) 



Here, £,„ is A^— dimensional vector according to the def- 
inition (|43|l . additionally each one of its components is 
a L dimensional vector due to its dependence on the 
Landau index Eq. (|40|) . In is the unit matrix and 
Dn — Diag {0,1,2, N — 1). Likewise, each of the 
components of the matrices Qrm In, and Dn are L x L 
matrices relative to the Landau contributions. In rela- 
tion with Eq. H28|l it was mentioned that the quantity 
g-zcrbfei does not explicitly depend on the longitudi- 
nal pseudomomentum component ki, this can be easily 
verified if the substitution c'„^ — e^^°'^''^ c,„ is incorpo- 
rated into Eq. (|^ . then using Eqs. (001 and (|^ . it is 
straightforward to verify that is indeed independent 
of fci. 

Eq. H45|l is one of our main results, and deserves to be 
emphasized. It is a generalization of the Harper equation 
that exactly includes the following effects: (1) an arbi- 
trary periodic potential, (2) an electric field orientated 
in a commensurable direction of the lattice, and (3) the 
coupling between different Landau levels. So far, no ap- 
proximations are involved, consequently it holds under 
most general conditions. In practice, only a finite number 
L of Landau levels can be included on the calculations, 
but a very good convergence can be obtained with a rea- 
sonable small selection for L. Previous known results are 
recovered if some approximations are enforced: («) if the 
electric field is switched off, its orientation is meaning- 
less, so the integers {mi, 1112) can be set to mi = 1 and 
7712 = 0, in this case Eq. 14511 reduces to the model pre- 
viously discussed by Petschel and Geisel^S. (n) li E — 
and additionally the coupling between different Landau 
bands is neglected and the potential is taken as the sum 
of cosines given in Eq. then the system reduces to a 
set of Harper equations, one for each Landau band. 

Let us return to the case that includes the electric 
field. Equation H45|) describes the system dynamics un- 
der most generals conditions; in spite of its complicated 



(44) 



structure the equation can be solved by the matrix con- 
tinued fraction methods. However, a considerable sim- 
plification is obtained if the electric field is aligned along 
the a;— lattice axis, i.e. (mi = l,m2 — 0). Henceforth, 
we consider that the electric field is orientated along the 
x— axis (mi = 1, m2 = 0), so 6 = a and additionally that 
the periodic potential takes the form given in (Q. The 
dimension N of the vectors Cm and matrices in Eqs. H43|l 
and (|44() reduces to = max{{r, s)(mi-f m2} = 1, hence 
the dimensions of Cm = Cm and Q reduce to include only 
the Landau indexes. Using equations H4U|I . (|41|l and H44f) 
'the matrices Q take the form 



= {Q+Y" = 
ttK 



{Q 

{Qmr = 



XttK 



2(1 + A) a 



2(1 + A) a2 

_ ^i{2'K a 'm-\-(Tbk 2 ) jjl-ii^ 



) + a 



(46) 



here the parameter K is a. measure of the strength of the 
coupling between Landau bands 



K 



ma^Uo (1 + A) 



(47) 



Taking into account these simplifications and that Dn 
cancels out, we find that Eq. I|45(l reduces to 



"C„i_i + QmCm + Q^C,n+l = (f + Ek2 + (jhEm) Cn 



(48) 




a 



FIG. 1: The energy spectrum inside the lowest Landau level 
as a function of the inverse magnetic flux a. The energy axis 

is rescaled to e = f / I^C/oe"'"'/^] . AU values of fci G [0, 2^/176] 

are included. The parameters selected are: m* = 0.067me, 
a = 100 nm, l7o = 0.5meV, E = 0.5l^/cm, and k-i = 7r/2. 



B. Numerical results 

The generalized Harper equation is given by a tridi- 
agonal infinite recurrence relation H48|) . If E is switched 
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0.4- 
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FIG. 2: The energy spectrum for the four lowest Landau 
bands is plotted as a function of magnetic flux a. The cou- 
pling strength K = 6 yields important Landau mixing. All 
values of fci £ [0,2n /qb] are included. The parameters se- 
lected are: m* = 0.067me, a = 100 nm, Uq = 1.2meV, 
E = 0.5V/cm, and ^2 = ■k/2. Seven Landau levels are in- 
cluded in the calculation in order to attain convergence. 



off, the equation becomes periodic with period p, in this 
case the equation can be recast as a finite pL x pL matrix 
that can be solved by a direct diagonalization. However, 
the introduction of E breaks the periodicity, and most 
of the methods used to solve the Harper equation break 
down. Eq. H48() was solved using an expansion of the as- 
sociated Green's function into matrix continued fractions 
(MCFs)21. The energy spectrum is determined detecting 
the change of sign in the Green's functions that appear 
in the vicinity of a pole. The density of states can be 
obtained from the expression 



1 



N {£) = -Tr G {£) = T-Im [Tr {£)] 



(49) 



were the discontinuity in Green's function is given by 
G{£) = G+{£) - G- {£), and the retarded and ad- 
vanced Green's functions are defined on the side limits 
G^ {£) = lim^^o G {£ ± le) . The numerical solution is 
obtained by truncating the iteration of the (MCFs) after 
the M— th term. The solution converges if M is large 
enough, in the calculations we observe that in order to 
obtain a convergence with a precision of one part in 10^, 
the cutoff can be selected as M fa 100. 

In our calculations we have used the effective mass 
m* = 0.067me typical for electrons in GaAs and a 
superlattice constant a = 100 nm*. The rescaled en- 
ergy spectrum £/ [/7oe~'^°'/^] is shown in Fig. Qfor the 
lowest Landau level, a weak modulation is considered 
(Uq = O.BuieV) so the ^ — v Landau mixing is negli- 
gible. The electric field intensity is i? = 0.5V/ cm cor- 
responding to a ratio of the electric to periodic energy 
p = eaE/Uo — 0.02. In the strong magnetic region 
a e [0, 1] the Hofstadter butterfly is clearly depicted. A 
distorted replica of the butterfly spectrum can be still ob- 
served in the region cr e [1, 2]. As the magnetic field loses 
intensity, the effect of the electric field becomes domi- 
nant, the butterfly is replaced by discrete levels sepa- 




FIG. 3: Energy density plot for the two lowest Landau levels 
as a function of the electric field intensity { p — eaE/Uo ) for 
a = 1/2. All values of fci £ [0, 2TT/qb] are included. 



rated in regular intervals, i.e. a Stark ladder. Although, 
p — eaE/Uo is small, there is a regime in which the elec- 
tric field dominates over the periodic contribution, we 
can trace down the origin of this effect to Eqs. H^^8(l and 
()46|l showing that the periodic contribution is modulated 
by a factor e~'^°'/^. Consequently, in addition to a small 
Landau mixing, the condition to preserve the butterfly 
spectrum can be stated as 



eaE <C Uoe 



-7rCT/2 



(50) 



This condition restricts the intensity of both the elec- 
tric and magnetic fields, these are estimated as (insert- 
ing units): E < Uo/{ea) ~ 100 {Uo/meV)V /cm and 
B > TTh/{ea?) ~ 0.6T. Fig. |21 shows the energy spec- 
trum as a function of a when a higher coupling strength 
[K = 6) produces Landau mixing. For strong mag- 
netic field the effect of the electric field is small and the 
spectrum is very similar to that previously obtained by 
Petschel and GeiselS^. Three regimes can be identified in 
the plot: (z) For strong magnetic field (small a) the Lan- 
dau bands are well separated and butterfly structures are 
clearly identified, (m) In the intermediate region (a 1) 
the periodic potential induces Landau level overlapping. 
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(Hi) as a increase very narrow levels appear as a result of 
the electric field. The development of quasi discreet level 
induced by the electric field is displayed in Fig. |21 here 
the energy spectrum as a function of E for the first two 
Landau levels is presented. For a — 1/2 every Landau 
level splits in two bands, these bands evolve into a series 
of quasi discreet levels as the intensity of E increases. 



= 1/2 



= 3/2 






p= 0.002 




p=0 



FIG. 4: Energy density plot for the two lowest Landau lev- 
els as a function of the transversal pseudomomentum k2. 
Three values of a = 1/2,3/2,5/2 are selected, so two bands 
appear when E = 0. Three values of E are considered: 
p = eaE/Uo = 0, 0.002, 0.01. 

In Fig. ^density plots of the density of states as a func- 
tion of the k2 pseudomomentum and the quasienergy £ 
are presented. The density of states was calculated from 
Green's function discontinuity obtained from the contin- 
ued fraction expansion. Three values oi a = 1/2, 3/2, 5/2 
are selected, consequently two bands appear when E = 0. 
In the superior panels it is observed that both as E in- 
creases or B decreases, the spectrum for a fixed value of 
^2 evolves into a set of narrow bands. When ^2 is varied 
the spectrum becomes almost continuous, except for the 
small gaps that open between the bands. 

Based on the previous results it follows that the EMB 
model preserves the band structure. For weak electric 
fields the bands are grouped forming the " butterfly spec- 
trum" ; whereas as the intensity of E increases the bands 
are replaced by a series of quasi discreet levels: a "mag- 
netic a Stark ladder" . 



C. Adiabatic approximation 

In this section we exhibit an alternative expression for 
the effective Schrodinger equation that governs the dy- 
namics of the system. Based on this formalism we shall 



find an approximated adiabatic solution, that throws fur- 
ther insight into the physical results obtained in the pre- 
vious section. Let us consider the generalized Harper 
equation l|48|l . taking into account that the matrix Qm is 
periodic in m with period p, we find convenient to define 
the unitary transformation 



[y,aHm/p] J^^^^^^j^j , (51) 



where [m/p] denotes the integer part of the number, and 
the delta enforces the condition / = m, mod[p\. Notice 
that whereas the index m in Cm runs over all the integers, 
the corresponding label of the state di {4>) takes the values 
? = 0, 1, 2, . . . ,p — 1. The new vector state di{(j)) satisfies 
the periodicity condition 



-) 

qa; 



(52) 



and the transformation matrices fulfill the following prop- 
erties 

oo 
m— — OO 

E / " ^-^^ (^) (0) = 5m,m' ■ (53) 



Applying this transformation to the generalized Harper 
equation H48|) yields 

Hm (fci -I- 0, fe) d (fci + 0, k2) = 

£ + Ek2-iE^d{ki + ^,k2), (54) 

where the Hamiltonean is reduced to a p x p block form 



HM{ki+<t>, ka) 



/ Qo Q+ 
Q- Qi Q+ 
Q- Q2 Q+ 

V Q+ 



Q- \ 





Qp-i I 

(55) 

here Q = Q (fci + 0, A:2), Qm = Qm (ki + 4>,k2) and 
= [ki +0,^2). Taking into account that each 
block in (|55|) is a L x L matrix determined by the Lan- 
dau indexes, we have obtained that the transformation 
(I51|l reduces the infinite dimensional representation of 
the Schrodinger equation to a pL x pL representation. 
The price is that a non-local derivative term respect to 
parameter (/> has been added. In the absence of electric 
field, Eq. H54|l represents a finite dimensional eigenvalue 
problem, that can be solved by direct diagonalization. 
We observe from Eqs. (|54|l and H55|l that appears to 
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be directly related the longitudinal quasimomentum fci, 
in fact if we redefine ki + (p ^ ki in Eq. (|54|l yields 

HM{ki,k2)d{ki,k2) = {£ + Ek2-iE^^ d{ki,k2). 

(56) 

This expression confirms the form of the Schrodinger 
equation (|30|1 previously discussed, in which the "dynam- 
ics" is determined by a differential equation with respect 
to the longitudinal quasimomentum. However we now 
have an explicit finite-dimensional matrix representation 
for the Hamiltonean. 

Let the "instantaneous" eigenstates of Hm be /i^"-* 
with energies A*^"^ , i.e. 



i7M/i(") (fci, A:2) = A(") {k^M) h^"^ (fci, fc2) 



(57) 



where a labels the band state of the non-perturbed prob- 
lem. A solution of Eq. (|55|l can readily be obtained in 
the adiabatic approximation as follows 

{ki,k2)^h'^"^ (fcl,fc2)x 

4^1 + 4 / ' d(l)A^"^ (0, fe) + Z7^"^ (fci, fc2) 



exp 



E 



E 



(58) 



where the Berry phase 7 (0) is determined from the sub- 
stitution of the previous expression in Eq. 



7^"^ (fci,fc2)= ^ r'd0/i(")t('/',fc2)|-/i^") (0,^2) 



(59) 

The energy eigenvalue is then determined by the period- 
icity condition (|52|) , hence the change of the phase of the 
wave function 158|1 must be an integral multiple of 2tt, 
consequently the spectrum in the adiabatic approxima- 
tion is given as 

nn /'27r/ga 

(fcs) =n£;aCT+ — / rffciA"(fci,fc2) 
27r Jo 



+ ^7M(27r/ga,fe), n = 0,±l, 
ztt 



(60) 



This result deserves some comments. The energy spec- 
trum in the presence of electric and magnetic field con- 
tains a series of discrete levels separated by 



AS = Eaa ■ 



hE 



(61) 



where we have used Eq. (|ll|l and restored units. These 
levels are similar to the Wannier levels that appear when 
an electric field is applied to an electron in a periodic 
potential, the energy separation being proportional to 
aE. In the present case, the band structure parallel to 
the electric field is replaced by a set of discrete steps, 
this "magnetic-Stark" ladder are characterized by a sep- 
aration proportional to the electric field intensity, but 



inversely proportional to both the lattice separation and 
the magnetic field. The existence of these levels can be 
explained by the following argument. In the presence 
of simultaneous E and B fields, the electron travels be- 
tween lattice points in a time t = aB/E. As long as the 
electron does not tunnel into another band, the motion 
appears as periodic, with frequency to = 27r/T, corre- 
sponding to a series of energy levels whose separation 
A£ — hid coincides with the result in Eq. H61() . This 
magnetic-Stark ladder is combined with the Hofstadter 
spectrum represented by the second term in Eq. H60|) 
and the contribution of the Berry phase, the competition 
between these three factors was discussed in the numer- 
ical analysis of the previous section. We notice that the 
Berry phase is written as the integral of the longitudinal 
component of the Berry connection Aa (fc) discussed by 
KohmotoS, (see Eq. 



V. HALL CONDUCTIVITY 

In order to calculate the electric current for a system of 
independent electrons let us consider the fermionic field^S 

i;it,x) = Y,J (1^) fc.,k {t, x) , (62) 

expanded in terms of the EMB states <Pa, that accord- 
ing to the Floquet Eq. are related with the solution 



of the Schrodinger as '^a 



ipa- The creation and 



annihilation operators bj, (fc) and ba (fc) fulfill the usual 
fermion anti-commutation relations (fc) , foj^, (fc')| — 
Sa,a'S (fc' — fe) , and zero for other anticommutators. 
The current density operator tjj''Jtp takes the form 

J=T. [4^ [^bi{k')b^{k){a',k'\U\a,k), 



J (27r)2 J (27r) 



where 



(63) 



(a',fc'|n|a,fc) = J d^xipl,^kit,x)n^^^k{t,x). (64) 

In order to evaluate the previous expression, we first con- 
sider the matrix element of x. Using the replacement 
ip (fc; t, x) — e^^'^u (fc; x) , that follows from the Bloch 
theorem Eq. H33|) . it can be readily demonstrated that 

(a', fc' \x\ a, fc) = —i^kS (fc' — fc) 5a' 

+ 5 {k' - k) A"' {k) , (65) 

with the definition 

(fc) - / d^X <, (t, X) l^kUc.k (t, x) . 

(66) 
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The subscript means that the spatial integration is re- 
stricted to a magnetic unit cell. We observe that the di- 
agonal element Aa (k) — A""'" (fc) is the induced Berry 
connectionSiSil. The derivation of Eq. (|65|l for Bloch 
states is well known (see for example^), the demonstra- 
tion is easily extended for the MB states^ and also for 
the EMB states^. Using the relation 11 = i[H,x], the 
effective Schrodinger equation l|5n|l and the result in 
one finds 

{a',k'\n\a,k) = 

[n„ (fc) S (fc)(l-<S„.,„)]5(fc'-fc), (67) 

where 

n„ (fc) = Vfe£„ (fc) - • VkAo, (fc) , 
„ (fc) = « (fc) - fa (fc) + tE ■ Vfc] (fc) . 

(68) 

Substituting (|67|l into l|63(l the electric current becomes 
J = 51 / I^*^^" (fc) - • Vfc Ac (fc)] 

+ * E / 70 (fc) (fc) [^"' (fc) - (fc)] " (fc) • 

a, a' 

(69) 

At very low temperature the electrons occupy all the lev- 
els that are below the Fermi energy. So the state of the 
system in the independent particle approximation can be 
written in the form 

i^)= n n^"(fc)io)' (70) 

where |0) is the vacuum, fp is the number of filled bands 
below the Fermi level and the state is normalized in such 
a way that the number of charge carriers is given as Af — 
vp S/{qb'^), where S is the area of the sample. For the 
state (|70ll the electric current becomes 

= E / I^'^^" (fc) - ^ • ^kAa (fc)] . (71) 

The first term on the right hand side represents the usual 
velocity group contribution. Additionally there is a novel 
contribution arising from gradient of the Berry connec- 
tion along the longitudinal electric field direction. 

In IIV CI we discussed the solution of the effective 
Schrodinger l|56() equation in terms of the adiabatic ap- 
proximation. However, if we consider the matrix element 
of Eq. (|56|l taken between the exact EMB states that 
solve Eq. (|56|l . we find that the energy eigenvalue fulfills 
the relation 

(fc) = A„ (fc) + £; • Aa (fci, fca) - Ek2, (72) 

where 

Aa (fc) = {Ua,k\HM\Ua,k)- (73) 



The eigenvalues in (|72|l and the previous results in H60() 
look similar; however the adiabatic approximation makes 
use of the approximate instantaneous solution H58|l in or- 
der to evaluate the quantities that appear in Eq. 16U|I . 
Instead the evaluation of (|72|l requires the exact solution 
of the Schrodinger equation (|56() . Substituting Eq. I|72() 
into the expression H71(l for the current, this reduces to 
the form 

J = 51 / (Vfc [E • A. (fc) - Ek2] - E ■ VkAa (fc)) 

=7r-EeT - —E x ) / — fc) , (74) 

a<VF 

the Berry curvature IIq, = e^^a is related to the poten- 
tial vector Aa by 

na (fc) = [Vfe X Aa (fc)]3 . (75) 

Eq. H74|l shows that the longitudinal magnetoresistance 
exactly cancels, regardless of the value of E. On the 
other hand, the Hall conductance can be read from Eq. 
()74|l and written (restoring units) in a compact form as 

g2 I 

CT// = — CTa , (Tq = - (1 - g?7„) , (76) 

n ^ — ' p 

a<UF 

where rja is the a subband contribution to the conduc- 
tance given as 

77„ = i /^[f7„(fc)]3 = i/ ^•A.(fc), (77) 

here CMBZ denotes the contour of the MBZ and the 
Stoke's theorem was used to write the second equality. 
We find it convenient to multiply and divide by q the 
contribution of rja that appears in Eq. H76|l . the reason 
is that the Brillouin zone ki € [0,2TT/qa], k2 G [0, 27r/a]: 
can be divided in q restricted Brillouin zones, RMBZ: 
ki G [0,2'K/qa], ^2 G [2Trj/qa,2Tr{j + l)/qa] with j = 
0, g — 1. The wave function satisfies periodicity con- 
ditions that connect the solutions in different RMBZ and 
thus the integrals in Eq. (|77|l can be evaluated in any of 
the RMBZ giving the same resulfe^. 

An expression similar to Eq. H74|l was originally dis- 
cussed by Thouless et al^ and Kohmotc^. However, their 
work is based on the Kubo linear response theory, conse- 
quently only terms linear in E appear in that case. Our 
result reduces to that of Thouless and Kohmoto if we 
drop the electric field in Eq. H56|) . in other words if we 
evaluate the conductivity in Eq. H77() using the zero order 
solution Eq. H57|l . As far as the gap condition is satisfied 
the argument that shows that the expression in Eq. H77(l 
is related to the first Chern number can be applied to 
prove that rja is quantized. However as the electric field 
is increased the high density of levels makes it practically 
impossible to resolve the conductance contributions for 
each miniband, as far as the original Hofstadter miniband 
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n = 2 



n = 3 




fi = 2 



(1=1 



(1 = 



FIG. 5: Berry curvature distribution on the MBZ for the 
first three Landau levels, a labels the Landau level fi and 
the corresponding subband n. The inverse magnetic flux is 
a — 1/3, so there are three subbands for each Landau level. 
The coupling strength is K — 0.2. 



is concerned this fact can be evaluated as a nonlinear 
contribution to ctq, that will lead to a break down of its 
quantization. Higher order corrections to the Hall con- 
ductivity can be calculated by a perturbative evaluation 
of Eq. H5t)|) . We observe that the term E{k2 — id/dki) in 
Eq. (|56|) can be considered as a perturbative potential; 
hence a series expansion in E can be worked out. The 
Hall conductance for the a— band is expanded as 



.(0) 



(78) 



Starting with the zero-order solution of Eq. H5t)|) 
Ua^ {t,x) = (t, a^ja'-^-'), the leading contribution to the 
Berry connection is simply given by 



(0) _ 



Vfc a 



.my 



(79) 



Then based on definition l|77|l and the condition H37() 
fulfilled by the wave function on the boundary of the 
MBZ, it follows that 77^"'' is the change of the wave- 
function's phase around the integration loop. The wave 
function has zeros in the interior of the MBZ thus 
Tja^ must be an integer. r]a^ is evaluated numerically. 
The contour integration around the MBZ in Eq. H77|) 
is numerically unstable; thus it is more convenient to 
evaluate the surface integral of the Berry curvature. 
iFrom Eqs. ^ and ^ is written as rji"' = 

I [(9fc^a(°) ISfcjQf^"^) — ((3fc2a'^°^ |9fe^a(°^)] , inserting a com- 
plete state "^2/3 \ inside the dot product and us- 
ing the identity 

(4°'-4°')(»'°'l£'5'»')-(«'"'||£|^""')> 

(80) 

that can be deduced from the partial derivative of Eq. 



(|57|l . the Berry curvature can be written in the form 

,(0) |H ^(0))(^(0) |H ,(0)^ ■ 

— C.c. 



dky 



where copa = ^f'' — ^a' ■ Eq. l|RT|l is well suited for 
numerical evaluations. Fig. shows the distribution of 
the Berry curvature plotted in the magnetic Brillouin 
zone when the inverse magnetic flux is selected as cr — 1/3 
and the coupling strength = 0.2. In this case, each 
Landau level (l is splitted in three subbands, results are 
presented for the first three Landau levels. It is verified 
that yyi^^ always yields an integer value; here a — (/i, n) 
labels the Landau level (i and the n internal miniband. In 



(81) 



and r^^^l = 1, 



this example we find: 77^*^^ = 1, 77^*^2 — 
regardless of /i. The numbers q y p are relatively prime, 
so there must be integer numbers u y v such that up + 
vq — 1. We can identify v = ria \ hence (1 — qvia^)/p 
is also a integer, but according to Eq. (|76|l this number 

is the Hall conductance of the subband a; thus Ua'^ is 
quantized in units of /h. As the coupling strength K 
increases some of the bands crosses (see Fig. ^ and the 

values of 77!°'' are exchanged. Fig. [7|shows the plot of fji"' 
for cr = 1/3 and K = 0.28. The Berry curvatures for the 
first Landau level are modified, now: = 1, rlf}^ = 1, 

and = —2. The profiles of fi^''' for small values of 
K coincide with those that appear in the semiclassical 
dynamics approach introduced by Chang and Niu'^^, our 
results include the effects produced by Landau mixing. 

Tabled shows the values of 77!° ■* for various selections 
of a when K is small. If the coupling between Landau 
levels is small, the following sum rules are satisfied 



0, 



(0) 



(82) 



here the sum includes all the subbands within a given 
Landau level. The last result guarantees that the Hall 
conductivity of a completely filled Landau level takes the 
value /h. However as the Fermi energy sweeps through 
a Landau level, the partial contributions of the sub- 
band, although integer multiples of e^/ft-, do not follow 
a monotonous behavior. If the Fermi energy lies in the 
71-minigap, the accumulated conductance C,a is defined as 
Cm," = Sj=i ■ "^^^ ri-minigap conductance satisfies 
the Diophantine equation n = Ag + pC^,"' (l-^l — in 
agreement with the results obtained by Thouless et al.^ . 

We now turn our attention to the higher order cor- 
rection to the Hall conductance. ^From the first order 
corrections to the wave function, the corresponding con- 



13 




FIG. 6: Energy density plot as a function the coupling 
strength K for the three minibands of the ji = 1 Landau 
level, a — 1/3. Notice that at 7^ = 0.28 the two upper bands 
crosses and the corresponding values of rj^^ are exchanged. 



n=l n = 2 n = 3 




tribution to the conductance is worked out as 



E 



dk 

CMBZ 27r 



«(o) 1111^(0)^^^(0) |||L|«(o) 



C.c. 



(83) 



As previously discussed, the wave function evaluated at 
opposite edges of the MBZ differs at most by a phase 
function fi (k) (see Eq. (|37|l 'l. this phase exactly can- 
cels in Eq. hence ^ = and there is no 0{E^) 
correction to the Hall current. 

The evaluation of rja requires the first and second or- 
der wave function contributions, then it leads in princi- 
ple to a lengthly expression, however the result simplifies 
if we notice that any term in which the wave function 
cancels will not contribute, because the corresponding 
contour integral around the MBZ will vanish. Then in- 
cluding only the relevant contributions the second order 
correction to the Berry connection can be worked out as 



(0) 



(84) 



(2) 

In the numerical evaluation of r]a we find convenient 
to scale the energy as e = [j/ge '^'^^^j : then as dis- 
cussed in IIVBI the spectrum for a given Landau band 
(for E — 0) lies in the region e E (—4,4). In this case 

Tja^ takes the form r/a'^ — 2Tr (elo/Uo)'^ e^"'fia \ where 
fja^ is dimensionless. The contour integration of l|84|) 
around the MBZ yields the second order correction to 
the Hall conductance. The expectation values of the par- 
tial derivatives of the Hamiltonean are numerically eval- 
uated in a similar fashion as those terms in Eq. (|81|l . 
As already mentioned the relevant part of A!"^^ is the 
phase, that integrated around the MBZ yields the lead- 
ing order contribution to an- Although the numerical 



FIG. 7: Berry curvature distribution on the MBZ for the 
first three Landau levels, a labels the Landau level and the 
corresponding subband. a = 1/3, so again there are three 
subbands for each Landau level but now K — 0.28 . Notice 
that the curvature distribution for the bands of the first Lan- 
dau level are modified as compared with those in Fig. |S] The 
reason is that K lies to the right of the point (see Fig. |SJ at 
which the second and third minibands crosses. 



evaluation of this phase is unstable, it was previously ex- 
plained how to determine its integral around the MBZ, 
some number are quoted in Table We can now ex- 
ploit a particular gauge invariance of the problem. Sup- 
pose that Ua^k satisfies the Schrodinger (|551) . then so does 
Ua,k — e^^'^'^^Ua^k, where g (fc) is a smooth function of ki 
and k2- A valid gauge transformation, must not be sin- 
gular, in such a way that its topological signature {rj^^^) 
remains untouched. Then, along a contour F, the phase 
of Ua,k can be choose in such a way that the Berry con- 
nection takes the form A!-^^ — 2imia''h{k) -f iVkp{k,), 
where h{k) is a selected function of k which integrates 
around the contour of the MBZ to 1. The term V^p (fc) 

does not modify the value of via^ , in particular it can be 
adjusted to have the same value in opposite sides of the 
MBZ. The gauge transformation cannot be implemented 
in all the magnetic Brillouin zone, it will be valid only if F 
encloses all the phase wave function singularities (zeroes 
of Mq); we are interested in the case in which F coincides 
with CMBZ. With all these considerations the second or- 
der contribution to the Hall conductance can be recast 
as 



= \ ' ^ dk- 

CMBZ 



dk 



(85) 



A suitable function h (k) can be selected as the gauge 
potential of a unit flux tube in fc— space: h (k) = 
^(fci, — A;2)/(fc? + fcf). In Table M we quote some of 

the values obtained for fja for the selected Landau sub- 
bands. It is observed that as the internal structure of 
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(7 


= 1/3 






= 1/5 




(J 


= 1/7 




n — 
H 


1 23 — 3 




/7 — 

y 


1 n — 5 






1 T> — 7 


n 






n 






n 




fjr 


3 


1 


0.416 


5 


1 


-0.018 


7 


1 


0.0239 














6 


1 


0.2760 








4 


1 


-4.152 


5 


1 


-110.2766 


2 


-2 


-0.328 


3 


-4 


17.836 


4 


-6 


117.3056 








2 


1 


-13.660 


3 


1 


-7.3303 














2 


1 


0.0011 


1 


1 


-0.087 


1 


1 


-0.006 


1 


1 


0.0003 



TABLE I: Zero and second order coefBcients 77^ , and ija 
{a = {fi,n)) for a — 1/3, 1/5 and 1/7. The miniband conduc- 
tivity aa is obtained from Eq. 17611 . The coupling between 
Landau levels is small, so the results are the same, regardless 
of the /i-Landau level. The sum rules (1821 and 1861 are easily 
verified. The coefficients are symmetric respect to cr = 1/2, 
the values for a ^ 1 — a are obtained with the replacement 
T]^'' — » —rja^ and 77!^' ffa^ respectively ; however Oa does 
not share this symmetry . 



band increase, the corrections to the conductance of the 
internal minibands increases abruptly. It can be easily 
proved from Eqs. (|84|l and (|85|l that the following sum 
rule applies 

E^M'n = 0- (86) 

n 

(2'] 

Consequently, even if each miniband presents cr^ oc 
c^/Uq B corrections to the Hall conductance, the contri- 
butions from a completely filled Landau level cancels ex- 
actly, and the Hall resistance quantization remains valid 
to order 0{E^). The nonlinear correction to an is then 
expected to be more important for the minibands. A 
critical electric field Ec can be defined as the the value 

(2) 

at which the second order contribution r]a becomes a 
fraction / of mia'^ , it is expected that the quantization of 
the Hall conductivity breaks down if the field becomes 
larger that Ec- According to Eq. l(75jl Ec is estimated as 



a = lOOnm and Uq = 0.5meV, the critical electric field is 
estimated as Ec ^ 215 V/cm for a = 1/3, Ec ^ 45 V/cm 
for (7 = 1/5 and Ec 23 V/ cm for cr = 1/7. Thus, the 
nonlinear effects are enhanced as the number of internal 
minibands increases. 



VI. SUMMARY 

In conclusion, we address the electric-magnetic prob- 
lem, and the calculation of the Hall conductance beyond 
the linear response approximation. We presented a thor- 
oughly discussion of the symmetries of the EMB prob- 
lems and of the construction of the wave function and 
effective secular equation. The dynamics of the system 
is governed by a finite difference equation that exactly 
includes the effects of: an arbitrary periodic potential, 
an electric field orientated in a commensurable direction 
of the lattice, and coupling between Landau levels. In 
addition to the broadening of the Landau levels induced 
by the periodic potential, the effect of the the electric 
field in the energy spectrum is to superimpose equally 
spaced discrete levels; in this "magnetic Stark ladder" 
the energy separation is an integer multiple of hE/aB. 
A closed expression for the Hall conductance, valid to all 
orders in E is obtained. The leading order contribution 
is quantized in units of e^//i. The first order correction 
exactly vanishes, while the second order correction shows 

a cr^'' (X /Uq B dependence. From the sum rule H86() 

(2) 

it follows that a)^' cancels for a completely filled Landau 
band. Hence the nonlinear correction to the conductance 
is expected to be of order 0{E'^) for a filled miniband, 
whereas for a complete Landau level the correction is 
expected to be of order 0{E^). We find that in order 
to preserve the self-similarity structure of the Butterfly 
spectrum, the electric field is restricted according to Eq. 
(|50|l . and additional restriction arises from the condition 
that the nonlinear correction to the Hall conductance re- 
mains small as given by Eq. (|87|l . 
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